Measuring the squared magnetisation of a molecule 
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■ It is proposed to shoot a molecule through a pipe of mesoscopic scale. The molecule should 

' effectifely carry two spins and antiferromagnetic chains should be embedded in the wall of the tube. 

If the momentum is detected before and behind the tube, a projection of the molecular spins B\ and 
<?2 takes place. If the molecule did not loose momentum inside the tube, it has a vanishing squared 
C"| ' total spin in the direction of its anisotropy axis /. Whereas [Z(<xi + f?2)] 2 = 1 if it did so. 
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I. INTRODUCTION 



It is one of the most broadcasted understandings from Quantum Mechanics, that a complete 
^ ■ measurement - e.g. the simultaneous, arbitrary well detection of the position and momentum of a 
^ ■ particle - is impossible. Thus one could be lead to the opinion, that quantum mechanics provides 
less principle possibilities than a classical model. That this is not true can be seen the most clean 
in the ideas of quantum computation 0,0]. To push these formal ideas into realization, one has to 
get control over the way the wavefunction gets processed. 

It has been shown that an important elementary step in this processing is an quantum mechanically 
^ \ incomplete measurement - but not in the sense above. It is easier to understand the distinction 
by means of discrete particle properties than with momentum and position. We use the working 
q ■ horse hilbert space of two spin half particles for it. The spin part of the wavefunction can be in a 
superposition of two spins up, | ft), two spins down, | f |), and the mixed versions | f |), | |f) jX 31 ] . 
A quantum mechanically complete measurement would now provide to the measurer one out of four 
possible results. The pointer of the detector shows e.g. 1 for | f f ), 2 for | ||) etc.. After that, the 
wavefunction is left in the detected state. It was complete, because the measurer could not have 
^sO . obtained more information about the particle. Determining the total magnetisation in z direction, 
thus 1,0 or —1, would e.g. be an quantum mechanically incomplete measurement. The difference 
is, that for the outcome 0, the wavefunction is still in an unknown superposition of | f f) and | | f). 
The measurement that was proven to be useful for quantum computation is another one 0, 0] • It is 
the Sf Qt = [/((?! + O2)] 2 ( or more generally the parity-) measurement along some axis I. In each of the 
two possible outcomes, the wavefunction is still left in a superposition. There are several proposals 
for this kind of measurement 0, 0, 0, H| . All of them have in common that the spins (which contain 
£3 ■ the quantum information) are localised in quantum dots. A macroscopic current can pass by, acts 
as detector 0] and projects the dots through coulomb interaction into the Sf ot = or S^ ot = 1 space. 
Our paper complements those proposals in the sense that the information carrier is moving and the 
detector is stationary. 



II. SETUP AND REQUIREMENTS 



The here proposed experimental setup is sketched in figure It consists of four distinctive parts. 
The magnetic oven provides a beam of molecules, that are each in an unknown superposition of 
their internal magnetic states | ff),| ff),| II), | ft)- We write them in the laboratory fixed basis 
in z-direction. They are assumed to leave the oven non entangled with any dynamical degree left 
in the oven. They are further in a product state with the centre of mass of the molecule. The 
splitting of their energies fl a is besides the kinetic energy of the molecule itself large compared to 
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all energy transfers in the further discussion. Once a while some molecule may pass the velocity 
selecting wheels and a state 4>R At k A {k) with mean momentum hk A is prepared at time t A at position 
Ra- It passes through the pipe, can excite the antiferromagnetic chains embedded in its wall and 
becomes ideally detected at time t B in a wavepacket 4>R B k B {k) characterised by k B ,R B and a B . 

<j>x{k) = (2no x ) {llA) exp[-al(k - k x f + i(k x - k)R x ] 

X = k x ,R x (1) 

The nature of the preparation and the detection determines (ta,b- If one uses e.g laser barriers to 
determine the time of flight , different light wavelengths can project the molecule into wavefunctions 
of different shape. 

We will see in the following, that the spin channels | fl), I IT) can n °t undergo inelastic scattering, 
whereas the channels | ft), | j j) can loose some of their kinetic energy and excite the antiferromagnet. 
An almost good setup for a parity measurement requires for the transmission amplitudes 

|I3I « ra forx=|TT>,|U>- ( 2 ) 
The observed momenta Tik B for the inelastic scattered molecules are expected to fulfil 

(k A -k B ) 2 > — + — • (3) 

Thus a momentum detection projects into the S^ ot = 1 space , if the momentum change is much more 
than %^Jl/a\ + l/a B and to S^ ot = 0, if it isn't. The word almost is referring to the projection to 
Sf ot = 1. Because the reacting channels can each excite different states in the antiferromagnet, there 
could in principle be situations in which the further processing of the molecule after the detection 
at Rb can not be explained by just looking at a superposition of the | ]]), \ [[) states. In the 
worst case, the spin degrees are entangled like | ]])\EX^) + | |J,)|.EX|) an d the two excitations 
inside the tube, \EXy±), decay each in independent processes to orthogonal bath states (e.g.' lost 
phonons'). The remaining molecule spin states behave then just like an incoherent mixture ??. But 
a quantitative discussion of that point is deferred to chapter HV1 and we first give the inequation (0) 
a microscopic foundation and elaborate how to design the experiment to fulfil it. 



III. EXPECTATION VALUES AND TRANSMISSION AMPLITUDES 



The expectation values to detect in the channel x the molecule with momentum hks under the 
condition, that the target stays in the groundstate or the n-th chain is left in the A-th excited state, 
are related to the transmission amplitudes by 
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FIG. 1: Sketch of the proposed squared magnetisation measurement. A molecule with spin quantisation axis / (along the large 
spaced dashed line) is emitted from the source at the left. Its spins <ti,c?2 are in an arbitrary superposition. It has to pass 
velocity selecting wheels and propagates in z direction (along the dashed arrow). Then a first laser beam detects it at time 
tA at position Ra- It passes a tube with antiferromagnetic chains in its wall. Some part of the molecular wavefunction leaves 
a trace in the tube by exciting spin waves and looses kinetic energy. Then there is at time ts a space-momentum detection 
at Rb- All molecules which lost sufficient momentum got projected into the subspace with [Z(<?i + S2)] 2 = 1, the others have 
[f(ffi + a 2 )f = 0. 

r )\ n create magnon excitations in the antiferromagnet, they are defined in equation (jB6j) . a k is a 
creation operator, it is related to the field operator of the centre of mass coordinate of the molecule 

\R)(R\ = J dkdk'e i{k ~ k ' )R a\a kl . (8) 

Without loss of generality can it be chosen to be bosonic because molecule molecule interference is 
of no relevance here. The time evolution operator 

U^\t) = exp{-it{H^ + V$)/K) (9) 
depends parametrically on the index of the quenched molecular spin state x. In the hamiltonian 

Ht ] = n x + Jdk E(k)a\a k + W^ll^n 

E(k) = h 2 k 2 /2m (10) 



are only the chain degrees included to which the molecule spins couple by the magnetic dipole dipole 
interaction 



V$ = T x J2x, n J dk dk< ala kn { n V t x t(k - k') 

{1,-1,0,0}. 
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Appendix B derives this expression and shows that W is the magnon energy gap at long wavelengths. 
It will be of some help to split the potential into a part that is symmetric and antisymmetric under 
space inversion V tot (q) = V s (q) + V a (q). Only one magnon excitations are considered, equivalent with 
the Born Approximation. One can make use of the complex variables [ik,£,\,n, which label coherent 
states on top of the molecule-magnon vacuum |0), 



\fj) = exp ^2fJL k a£j |0), |0 = exp (j^^nliA |0). 



To perform the expansion to first order in the reaction probability of a single chain 

9° k At^) = (» k (t)iil<(0))o = 5 ktk ,G(k,t)F(t) 



dq dq l drV t X l n (q - q l 



X(/i fc (0^(T)/i*(T)//*(0)eA,n(t)elnW)o 



(12) 



(13) 



(14) 



one uses the Wick theorem to express the path integrals [111 



(••}( 



D\nn*££*\...exp 



Jdt \Y,F-\t)^ n (t)il n (t) 



dkG-\k,t){i k (t)n* k (t) 



(15) 



in terms of G(t, k) = exp(—ihk 2 t/2m)0(t), F(t) = exp(—iWt/h)Q(t). 0(..) is the unit step function. 
The elastic contribution describes just the projection onto the free propagating wavepacket 



m = x(tB-t A )- 1/A exp 
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-(k A - k B f — + — 

V °A a B, 

(K el - V el (t B - t A )f " 
<%! + a 2 B ) X (t B - t A ) _ 



V - R d V _ h k ^A + k B <J 2 B 

l<-el — H-B — -ft A Ve 



m a\ + a B 

x (t) = 1 + t 2 /t 2 dlsp t dlsp = 2V2m(a 2 A + a%)/h (16) 

and leads for all spin channels to a detected momentum that differs only through the detection and 
preparation process itself from the initialised momentum. The detection times are in the following 
assumed to be always much smaller than the generalised dispersion time tdi sp , so, that the molecule 
is at detection still in the regime of corpuscular behaviour. 



The propagation kernel for the inelastic contribution depends on the lower(upper) momentum cut 

off ^Q<(>) 
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-i5 k , k <V s x > n (0)G(k,t)(Q > -Q < )T x J dt l F{t-t l )/h 
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-iE(k')t/h 



-i(E(k)+W)t/h 



(17) 



E(k l ) - E(k) - W 

If one inserts this into equation (J7J), the first part gives rise to the same k B distribution as in equation 
Ijlfij) . Thus such a momentum detection wouldn't reveal the excitations inside the tube. We shouldn't 
worry about the fact that it even diverges if one relaxes the cut offs. Including an infinite momentum 
spectrum is for the wavepacket and the magnetic field anyhow unphysical. We simplify for large 
enough detection times, h < Min {E{Q < ) + W - E(Q > ), W) (t B - t A ), 



V s X ' n (0)(Q> ~ Q< 
W 



E 

A, 71 

+ J dkdk l (j) B (k)V t x f(k - k l )exp [-i(E(k) + E(k l ) + W)(t B - t A )/{2%)} 
2 sin [(E(k') - E{k) - W)(t B - t A )/2h] 



E(k)-E(k l ) - W 



^(fc')e 



X = | TT), | II) 



by substituting gm ((*a~^)--) with a proper normalised step function between —it/{t B —tA) and ir/ {t B — 
t A )- Time dependent phases, ipi n , appear. But they are independent of A and n and will through 
the following discussion become irrelevant. We are left with the task to perform the integral 

exp[-iE(k)(t B -t A )/h) 



dkv£T{k-Jk* + Q 2 )Mk) : 



k 2 + Q 2 ) 



(19) 



Vk 2 + Q' 2 

using the notation Q 2 = 2mW/h 2 . All the Besselfunctions in equations ()B11IB11|) behave for large 
momenta like exp(—\qd\). Thus if one wants to take full advantage of the potential, one obtains an 
upper limit for the tube diameter 
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(20) 



which can be ignored for k 2 B > \-^a B )- 2 • How to calculate the integral depends essentially on the 

ratio of a a and a B . Thus the approximate gaussian contribution around k = \jk\ — Q 2 and k = k B 
gives in the limiting cases 
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with the abbreviations 



Qin = k B - ^k A - g 2 v in = ^k A - g 2 

Kin = R b -Ra(1- Q 2 /k 2 A ) 1 ' 2 (tab = a A (l- Q 2 /k 2 A ) 1/2 



GST) 



for a A ^> a B and 



Qin = [k A - ^i + g 2 ] (i + g 2 AI) 1/2 °ab = o A (i + g 2 Air 1/2 
^ m = i? B - r a (i + g 2 /4)" 1/2 v in = ik B pi") 

in the opposite limit. If one ignores, that in equation (|21jl the sums over A, n are inside the |..| 
signs, the most probable momentum to detect is in both cases for these inelastic scattering paths 
h(k A — g 2 )" 1 / 2 . But what is now the effect of taking the sum first ? Important is the sum over 



the z-positions of the different spin rings r n .One gets \Tf n \ oc Y.n ex P(~^(k A — Jk A — Q 2 )r n ). If 
there are many spin rings these phases would usually average out each other and we don't find 
inelastic scattering at all. For this reason we have to design the distances of the rings so that 

r n = 27rJ\f(n)/(k A — \Jk\ — Q 2 ) , where A/"(..) gives only integer numbers. This resembles the idea 
of multilayered mirrors which achieve very high reflection for a certain wavelength of light. Thus if 
one is able to design the positions of the rings in the prescribed way, then one can improve as well 
the accuracy of the momentum transfer. 

Equations (|18I19J) offer two strategies two diminish the inelastic contributions that would lead 
to measurements of k B ~ k A , t B ~ t A + TZ^V^ 1 . If one is able to keep the spin quench axis of 
the molecule along its flight direction (J3 = ir/2 in equations (jBll)Blljl ). the symmetric part of the 
interaction V s vanishes. We insert this condition into the requirement (J2J) at the optimal measurement 
values. We can further assume that the tube is much shorter than its distances to the locations of 
preparation and detection and obtain our main result 



1 <ti 4/i 0/ u 



2 MVSNJ 7 {W) 



d 2 v A h 



o A a B 



[1 - W/E(k A ))a 2 A + a 2 B 



as an estimate for the minimal number of spin rings M. The function 

F(W) = exp -arctanh (- 




1/2 

(21) 



(22) 



goes monotonously from 1 at W/SJ — >• oo to oo at W = 0. It may be tempting just to decrease W, 
so that M doesn't have to be very large. But the validity of the Born approximation requires that 
the detection probability of a single chain is small. The opposite of condition (}2~T]) has to be true for 
M = l,thus it gives a lower limit for W as well. 

In the realistic situation J = W, S = 1/2, and a B 3> & A , these requirements can be written as 



vumit < v A < Mvumit with v limit = 108 J 



va = 439,/^- (23) 



e(2-e)m 

, where we introduced the relative momentum loss e. The initial speed v A and vu m i t are given in 
m/s. The mass is given in atomic mass units (A.M.U.) and W in its typical units of meV. p is 
the distance between two spins in a chain and, like the other lengthscales d, a A , as, should it be 
inserted in A. Figure El shows possible v A ,e combinations for a given molecule mass. 
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m [ A.M.U. ] 



FIG. 2: Solutions for equation 1231 are shown for various relative momentum loss e. The dashed line gives the upper limit 
Mvumit for the initial velocity va, the lower limit is far below the lowest shown line. We have chosen lmeV for the spin wave 
gap and a tube with 10000 rings, take in each ring a nearest neighbour distance of p — 3A a tube diameter of lOOA - thus each 
ring contains about 105 coupled spins. The relative momentum loss is taken as e = 1Q X where the upper line corresponds to 
e = 10~' ! and the lowest one to e = 0.1, x changes from curve to curve by 0.1. The molecule mass is given in atomic mass units 
(A.M.U.). 



The second strategy to diminish the unwanted contribution requires that the molecule rota- 
tion is in a chaotic regime. This means one is allowed to ignore the actual angle trajectories 
and to average over their distribution. We assume a homogeneous one and just substitute 
\sin(f3)\, \cos(/3) \ — > 2/ir etc. . The ratio of unwanted to wanted contribution is then in the order of 
(Q> ~ Q<)^ v a/{^W). If we insert lA -1 as estimate for the momentum cut off (about 10 times the 
inverse size of a buckyball), the ratio becomes small if 

97 W > v A (24) 

is fulfilled - provided one uses the same units as in equation ()23|). Equation (J2"T|) has of course to be 
fulfilled simultaneously, but figure El exemplifies that this is possible in a realistic parameter regime. 
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IV. JUST AN ALMOST GOOD MEASUREMENT ? 



This section should give more quantitative arguments to the discussion in the context of equation 
(J2J). We elaborate how an almost good measurement can become a really good measurement. 

We have to understand more about the state in which a spin chain is left after | | J, J,) 
passed by. It is much more difficult to calculate these states than just obtaining the transmission 
amplitudes in the last chapter, because one has to follow the full quantum dynamics. We show now 
that it is for the current purpose not necessary. 

We restrict ourself to the spin ring at r n = and use the notation |A) = ^alGS) for the excitable 
magnon states. In the same spirit as we were using the Born Approximation, we can project out 
of the Hamiltonian Hq + V£ ip the states with more than one magnon. The effective Hamiltonian is 
then in bra-ket notation 

Heff = J2E(k)\k,GS)(k,GS\ 
k 

+ Y, S kM (E{k) + W)\k, X)(k, X\ ±V t x Q f{k - k')\k, X)(k' ,GS\ + h.c (25) 

k,k',X 

where the '+' corresponds to the | ~f )- an d the '-' to the | ||)-channel. The evolutions equal under 
the transformation |A) — > — | A) , A = 1,..,4. An initialised state (fi^\ Tt) + /ijjj ID^CS) evolves 
then into (1 - \Z\ 2 )^ (/i TT | ft) - (JL U \ U» \ex) + Z (/x TT | TT) + A* U | ID) \GS). Where \ex) is some 
excited state of the spin ring and \Z\ < 1. If this state decays independently of the molecular state, 
it provides just an arbitrary phase Vt and in the final state 



^(z + ^i-\z\^ex P (in)))\n) 

\GS) (26) 



+ /i U \Z-Tjl-\Z\*exp(in)))\ll) 

the spin ring factors out. Thus in the limit \Z\ —>■ 1, the relative strength of the coefficients is 
not even changed by the scattering. The crucial aspect in the argument is that the decay happens 
independently of the molecular states. Having a longer relaxation time than the typical time which 
the molecule stays inside the influence of a single chain should be a sufficient condition for its validity. 



V. SUMMARY 



We have proposed a new type of experiment to measure the squared total magnetisation of a 
molecule in the direction of its spin anisotropy axis. We have shown in which parameter regimes 
one can realize it. The experimental side meets several challenges. The first key problem is to find 
an appropriate molecule. It should have the property, that it carries the basis to create the Bell 
states by two spins and that the individual spin multiplets have big energy splittings. It is further 
necessary that one is able to build a 'magnetic oven' - a source of these molecules, where they get 
excited into superpositions of these states. The second big question is the actual realization of the 
tube. It is a nontrivial goal to fabricate such a mesoscopic object especially with small diameter but 
still many individual chains. But we believe it is a way to go, because improved technology can help 
to make use of the peculiarities in Quantum Mechanics. 
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APPENDIX A: SPIN WAVE APPROXIMATION 

The spin chain is described by the antiferromagnetic heisenberg model 

H AF =jj2S A ( V )S B ( V ) + S A ((p+l/N)S B ( V )+SA(S z A ( V )-S z B ( V )) (Al) 

tp=l/N,2/N,.. 

, where (p labels the unit cells, A, B stands for the sub-lattices and A mimics e.g. single ion anisotropy 
and stabilises the Neel state. Through the transformation 

SaW) = a W ~ S > S z M = S-b% (A2) 

S\{ V )=a^2S-aU^ S+(<p) = y/2S^^ p b v (A3) 

S A ((p) = \j~2S - ajpa^ , S^(ip) = b^^S - bjpb^ (A4) 

(N-l)/N 

a v = jV"3 Yl e i2nk{lp +m) (cosh{e k )a k + sinh{Q)^ k ) (A5) 

k=0,l/N,.. 
(N-l)/N 

b v = N~* e ~ i2 " kLp (cosh(e k )P k + sinh(Q)al) (A6) 

fc=0,l/JV,.. 

, . ^ , cos(kn) 

tanh(2Q k ) = ^ - K ^ (A7) 

Referenced to Holstein and Primakoff, is a mean field hamiltonian 

H AF = -2JNS(S + 1) - JN(2S + 1)A 

+ Y + PhPk + l) + 0(S°) (A8) 



Uk = JS^{2 + A/S) 2 - Acos 2 (kn) (A9) 

obtained, which contains the decoupled magnons. It corresponds formally to a large S expansion 
around the Neel Groundstate. 

APPENDIX B: MAGNETIC DIPOLE INTERACTION 

The magnetic dipole interaction between the spins of the molecule and the chain (in the plane 
with z = 0) should now be expressed in the same approximation. 

V fffVnnT ^ [ifoOft] [R(<p)SM] 

R(<p) = {dcos{2ir<p/N),dsin{2ir<p/N),R) (B2) 

contains the magnetic constant of the vacuum /i , the electron spin moment fi e , the molecular z 
component, its spin degrees <7j. x labels the A, B sub-lattices, d is the tube diameter. The molecular 



10 



— * 

spins are quenched along the axis in direction / = (cos(P)cos(a), cos(P)sin(a), sin{f3)). Thus one 
obtains for the effective interaction (x\Vdi P (R)\x') = S XtX 'V^ ip (R) 







V1S } (R) = -Vjff (R) = V asym + V sym 



(B3) 
(B4) 



V, 



asym 



N' 



3dR^ ( sin{f3 y os{2 V )S *(y) + sin{(3)sin{2it^)Sl{y) 

\R{v)\ b x,^ v 

+cos(P)cos(a + 2*^)S'{<p) 



V. 



sym 



\R{<p)\* 



J2 (cos(f3)cos(a)S*(ip) + cos(/3)sin(a)S»(<p) + sin((3)S^<p) 
-3d 2 non 2 e 



+ 



\R(<P)\ l 



X,f 



( cos(l3)cos ( a + 2ir^- ) cos(ip) 



N 



xS*(<p) + cos((3)cos(a + 2itj-)sin{<p)S y x {ip) + sin(/3)Sfa) 



(B5) 



We have split the parts that are symmetric and antisymmetric under R — > —R. If one expresses the 
spins by the magnon operators one finds, that the interaction couples only to four degenerate modes 



7 = { a l/N, Pl/N, ai-i/jv, Pi-i/n} 



(B6) 



with the energy of the anisotropy gap W = 2JAy 1 + A/ S + 0(1/N) 

3RdfiQfi 2 



V B 



'SN 



asym 



sin(/3) 



(R 2 +d 2 ) 



5/2 



cosh(6 )(jt - 7 |) + sinh(6 )(jt + 7s 

+h.c. +0(1/N) 



(B7) 
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Vsym — Hof^e] 



3d 2 



>SN 



cos(j3) (cosh(9 ) + sinh(9 )) 



(R 2 + d 2 f/ 2 2{R 2 + d 2 f/ 2 V 6 71 + 6 72 73 + 74 

+h.c. +<D(1/N) 

It is desirable to express now all in the second quantised form 

V1S } = E / dkdtala k ^l n [V s x ' n (k - V) + V a x ' n (k - fc')] + /i.e. 

A.n 



:b8] 



(B9) 



, in which (A, n) labels the A-th of the four modes in the chain with z coordinate r n . The fourier 
transform lead to the modified Besselfunctions of second kind B\..] 



V a L > n {q) = ie-^^ e 



>SN 



sin{p)cosh(9 Q )q 2 B[l, \qd\] 



V'^{q) = V^ n {q) = tanh(9 )V a ^(q) and V a ^(q) = -V^{q) 
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<SN 



V s 1 > n (q) = V s 2 > n (q)* = e- ir ^^ 2 eX j^-cos{f5){cosh{e Q ) + sinh{e Q )) 



2e~ 



B[l, \qd\] - e* a q 2 B[2, 



'SN 



V^{q) = -V s 4 > n (q)* = e-^^ 2 e ^cos((3)(cosh(6 ) + sinh(6 )) 

xe~ ia q 2 B[2, 



(Bll) 
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